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1. Atmospheric Heat Engines on Earth and Mars. A terrestrial 
micrometeorologist is constantly aware of the intimate connection between 
the surface fluxes of H 2 0 and of sensible heat. This connection is 
ubiquitous evidence of the Earth's atmospheric heat engine in action. Its 
character depends, inter alia, on the tightness of the linkage between the 
energy cycle and the cycle of H 2 0. H 2 0 is an important component of the 
working fluid because the temperature and pressure at the Earth's surface 
are such that H 2 0 is always present there in at least one condensed phase 
as well as vapor; and because the surface air's latent heat content is 
relatively large compared to its sensible heat capacity, with consequent 
correlation between energy transport and H 2 0 transport. 

The same micrometeorologist finds however, that the latent heat capacity of 
surface air on Mars is only about one ten- thousandth that on Earth. 
H 2 0-based latent heat fluxes represent only a trivial fraction of total 
energy fluxes on Mars: the linkage between H 2 0 and energy fluxes is 
extremely weak [1]. Unlike on Earth, H 2 0 on Mars is an ineffective 
component of the working fluid. C0 2 is the dominant component. 

Can we make these considerations quantitative? Equation [1] defines 0, 
which we propose as an index of the relative importance of latent heat 
transport by a particular component under particular surface conditions : 

0 - ^ [ 1 ] 
pc p 1 

Here p is the density and Cp the specific heat at constant pressure of the 
surface air, p v the surface vapor density of the component, and L its 
latent heat of phase change (evaporation or sublimation) . 

The quantity 9, with the dimensions of temperature, is the excess of 
surface equivalent temperature (2) for the particular component, over the 
surface air temperature. We may call it the equivalent temperature excess. 
A tight linkage between energy flux and the mass flux of a component will 
give large values of 9 under given surface conditions; and conversely. 

Figure 1 graphs evaluations of 0 for H 2 0 on Earth and for H 2 0 and C0 2 on 
Mars. On Earth 9 decreases from 93 K for a surface dewpoint of 310 K to 
0.1 K for a frostpoint of 225 K. On Mars 9 for H 2 0 decreases from 6.2 K 
for a frostpoint of 220 K to 10* 5 K for one of 150 K. Contrast the very 
large value for C0 2 on Mars, 9 - 692 K. The latitudinal distribution of 9 
on the two planets signalizes vividly their different meteorology. Earth's 
atmospheric heat engine is most effective and active in the tropics and 
least so at the poles. On the other hand, the Martian heat engine is most 
effective over the polar C0 2 caps, with 9 there more than 10 times that of 
tropical Earth; but it is ineffective elsewhere. 
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Figure 1. 

Variation of 0 with surface condensation 
point for H 2 0 on Earth and for C0 2 and 
H 2 0 on Mars. 


2. Diurnal and Annual Cycles of H 2 0 In the Martian Regollth. The looseness 
of the atmospheric linkage of energy and H 2 0 surface fluxes enables 
analysis of the equilibrium and flux of H 2 0 in the Martian regolith with 
the atmosphere entering only through surface boundary conditions on 
temperature, T, and H 2 0 vapor density. The regolith is so cold that we 
need consider only the processes of sublimation, vapor diffusion, and 
thermal conduction. In such systems ice behaves as a "wetting" capillary 
substance. Many of the concepts of terrestrial soil-water physics (e.g. 
3,4) with appropriate modification, carry over. Details are given in 
(5,6), where this approach was used to investigate model diurnal and annual 
H 2 0 cycles. The vector flux density of H 2 0, q, is described by the equation 


q - -DflVfl - D X VT [2] 

where the coefficients D#, D x are strongly nonlinear functions of both T 
and of the regolith volumetric ice content 0. 


D x increases very strongly (approximately exponentially) with T, so that an 
harmonic temperature wave imposed at the surface produces a net downward 
flux due to the temperature field. For equilibrium over the year this must 
be balanced by a net upward flux due to 0 increasing with depth. These 
considerations yield the equilibrium ice content profile. In principle, 
this profile is perturbed by diurnal and annual oscillations in surface 
boundary layers. These perturbations, in fact, prove to be trivial. 

The earlier work produced ice tables (regolith saturated with ice) at 
depths below the surface which were very sensitive to the surface humidity 
cycle. In the light of later information, the humidities used then were 
unrealistically high. The present calculations are for conditions 
approximating the mid- latitudes . They are: mean surface temperature 200 K 
with diurnal and annual semi -amplitudes of 30 K; and daytime surface 
frostpoint (DSF) from 190 to 200 K. The solutions give 0 in the range 
0.02 - 0.03, with a slow increase with depth down to about 12 m. There the 
areothermal heat flux becomes important, reversing the gradient of 0. See 
Table 1. 
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Table 1. Equilibrium Ice Content Profiles 


Daytime Surface Frostpoint (K) 

190 

195 

200 

Depth (m) 


Ice Content 

6 

0 

0.02497 

0.02569 

0.02648 

0.1 

0.02630 

0.02719 

0.02819 

1 

0.02803 

0.02924 

0.03090 

2 

0.02883 

0.03034 

0.03246 

6 

0.03074 

0.03319 

0.03994 

12 

0.03122 

0.03455 

0.04209 


Figure 2 shows the annual cycles of H 2 0 flux into and from the regolith for 
DSF's of 190, 195, and 200 K. Figure 3 shows, for a DSF of 195 K, the 
diurnal H 2 0 cycles for 12 representative sols. The notable point about 

these results is the very small value of the H 2 0 fluxes. Even for the 
largest DSF, 200 K, the total annual variation of condensed H 2 0 in the 
regolith is only 120 x 10* 6 m. This corresponds to variations of 0 of only 
+ 0.00006 distributed over the top 1 m of the regolith. The largest 
diurnal change in ice content (DSF - 200 K, sol 167) is 1.2 X 10* 6 m, 
corresponding to a change in 6 of only 0.00003 distributed over the top 
0.04 m of the regolith. 
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Equinox Solstice Equinox Solstice 



Figure 2. 

Annual cycles of H 2 0 exchange with 
the Martian regolith. Illustrative 
examples with DSF - 190, 195, 200 K. 
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Figure 3. 

Diurnal cycles of H 2 0 exchange with 
the Martian regolith. Illustrative 
example with DSF - 195 K. Note 
that 24 hars (Martian hours) - 1 
sol. Numerals on curves show sol 
number, with sol 0 at vernal 
equinox . 
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3. Influence of Areothermal Heat Field on Regollth H 2 0. We have already 
noted that the areothermal heat field interacts with the damped surface 
temperature waves in the unsaturated regolith leading to a maximum value of 
6 at about 12 m. Below this 6 decreases slowly to negligible values in 
depth . 

Interest attaches to the possibility of liquid water in the regolith at 
great depths. The depth to the 273 K isotherm, 

z 273 - £ (273-T 0 ) , [3] 

with T 0 the mean annual surface temperature (at z - 0) , X the mean thermal 
conductivity in 0 < z < z 2 73 , and Q the areothermal heat flux density. We 
adopt Q - 0.025 W m" 2 and take X - 0.728 W m* 1 K" 1 , consistent with the 
profiles of 6 of the examples discussed above. We thus find, for 
T 0 - 200 K, the value z 273 - 2126 m. Typically, the possibility of liquid 
H 2 0 in the regolith arises only at depths of the order of kilometers. 

If the regolith H 2 0 is in equilibrium with the present Martian atmosphere, 
its H 2 0 content below the 273 K isotherm seems likely to be so small that 
the H 2 0 may exist primarily in an adsorbed phase. 

Matters would be different, of course, if a mass of fossil ice-saturated 
regolith were connected to the 273 K isotherm. In the unlikely 
circumstance where the region below 273 K was undrained, one would then 
have water- saturated regolith underlying the ice- saturated mass. 
Underdrainage seems the more likely possibility. In this case the final 
profile of H 2 0 beneath the ice-saturated mass would be established by H 2 0 
vapor equilibration in the areothermal heat field. 
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Figure 4. 

Case of ice -saturated regolith 
with ice extending to 273 K 
isotherm. Volumetric ice and 
water content profiles above 
and below isotherm. 


Figure 4 shows the profile of ice and water content for an illustrative 
example. For a regolith with underdrainage, we see that, even granted the 
existence of fossil ice -saturation above 273 K, liquid H 2 0 can occur only 
at a depth of order of kilometers , in minor concentrations . 
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